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NOTi^TIONS 


G-eneral. Sjmbols 

: Gross -sectional area 

: Oross-sectional area at "base 
: I'londimensional area = A°/a° 

: Area of water moving with the dam 
: Nondiniensional angular velocity = ijr^ L/k^ 

: T/idfch of water 

: Fundamental critical damping = 2 a*^ f 

j 

: Damping constant 

C 2 : Families of physical characteristics 
; Pressure coefficient 
: Modulus of Elasticity 

S I 

: iTondimensional forcing function = 

; Shear modulus 
: Acceleration due to gravity 

: Pleight of water in reservoir basin 
: Second moment of area about y-axis 



Second moment of area atout 7-axxs at base 

Nondimensional second moment of area ahout 
y-asix = l/l^ 

-R-otary inertia about y-axis 

Robary inertia at base about y-axis 

bondimensional rotary inertia about y-axis 

= j/j„ 

Sliape factor 

■velocity of flexural -ware = 

Velocity of siiear wave = p° 

Height of the highest uionolith of the dam 
Bending Moment 

Fondiraensional bending moment = ML/EI^ 

Radius of gyration about y-axis 
Time 

Shear force 

2 

ITondimensional snear force = Yh /El^ 
Fondimensional linear velocity = 
fisplaceqient in z-direction 

Rondimensional displacement in z-direction = \v'/l 



z 


; Distance from base 


X 

• 

« 

londimensional distance from bai 

^tt 

4 

dround acceleration 

2. 

Dreek 

SvTibols 

a 

• 

= 1.0 + aV P° a'" 

P 

« 
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2 p® fX 
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= P/C 2 


• 
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A 1^ 3 

(kgAp 
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• 

■^0 

X 
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A 

• 

• 

Dauiping ratio 


/E : 

D^ndimensional bending stress 

•" 7 ^ 

« 

n mdimens lonal time = kp/l 

0 

• 

# 

Angle of taper with, x-axis 

it 

m 

m 

Bending slope 



g A„ l5 

- 7 /Og 

randamental frequency 
Damped frequency 
i^afaral frequency 

Mass per unit volume of learn 
Mass per unit volume of water 
Jump 


3, Mathemabical Symbols 


C J = Matrix 



Row vector 
Ooluuin Vector 
Determinant 


: Transpose of coluron vector 



ABSTRACT 


The transient stresses developed due to earbhcLaake 
plaj a ma3or role in the failure of structure due to dyna~ 
mic loading. In this study integration along bhe character- 
istics of the resulting eq_uations is used in numerical 
solution for the transient response of dam-like sbructure. 
The method is compared with known closed form solution and 
IS found to be in good agre^jment. 

The drm-like structure is treated firstly without 
water in reservoir basin and secondly with impounding water. 
The effect oC water on the dam is taken into accounb by 
virtual mass addition. It is found that Lhe maximum tran- 
sient bending stresses developed in the dam a.re higher than 
those obtained by the usual mode superposibion analysis. 

This aaalysis is applied for the dynamic response 
of the highesb monolith of the Koyna Dam subjected to 
December 11bh (1967) earthq_uake. The method is equally 
applicable to bhe dynamic analysis of stacks subjected to 


earth motion 



GHi\PT6‘R I 


INTRODUCTORY RDivIABXS 


1.1 OTRO DUCT ION 

Ari ear'cliqaake is a ground vibration phenomenon, 

T>xe ground vibration associated with an earthquake is of 
random type i.e. the peateof acceleration, both positive 
and negative, have various amplitudes and occur at various 
time intervaJ-s without any specific pattern. Structures 
situated on the ground are forced bo respond to these 
motions and therefore set into vibration. An eorbhnuoke 
produces only ground vibration, therefore exerts forces 
only at the base of a structure. The forces induced in 
the structure it a higher level are developed as a result 
of dynamic response o" the structuie. During an e-irtlrnake 
the .ground motion occurs in all directions horizontal and 
vertical. However, horizontal components are more important 
than vertical ones since the structure is usually strong 
enough in the vertical direction. 

In general, three techniques are a,vailable for the 
theoretical auialysis of transient stresses in elastic bodies 
caused by impact loadings such as ea.'thquake. These uiethods 
are; 



(i) Th.e mebiiod of mode superposition 
(li) The Laplace transform technique 
(iii) The method of characteristics. 

The method of mode superposition is often unsatis- 
factory since it utilises the normal mode oC vibration as 
predicted by the elementary theory and many modes are 
ofbeia required for rapid convergence of results. Also this 
technique tends to round off or average out the sharp peaks 
in the response. 

The L place transform technique gives a solution 
of closed foro type, hence it always demands a closed 
form representation of input motion durin , the total 
interval of disturba3ice. Moreover, this technique has its 
inher^ent inversion difficulties and frequently numerical 
method of integration has to be employed to evaloiate 
inverse iotegral of the transform. 

The a,dvantages oT the method of characteristics, 
employed in this study, are that, firstly, the problem can 
be formulated continuous one i.e, no discxetiza,tion 
0 L physical rnodel is needed. Secondly, it gives a simple 
description of wavefronts and finally, well known technique 
of numerical integration can be effectively used. 



1.2 HISTORIGAl BAC.G&ROOTD 

The method of characteristics was first used hy 
Leonard and Budiansky jl] to solve a strictly mathematical 
transient problem of a beam i.e, when the governing eq_ua- 
tions and initial conditions are homogeneous but the boun- 
dary conditions are nonhomogeneous, Hov*fever, solutions 
were obtained for a particular type of material for which 
both shear and flexural wave have the same velocities of 
propagation. Subseq^uently Jhasman \2 [ obtained the relation 
between the various stresses and displacement derivatives 
at the wavefronts due to step inputs for circular sheets 
and plates following the same general principle of Leonard 
and Budiansky, but did not solve the characteristic equations 
behind the v/ave front. Later Plass \ 3 \ improved the problem 
of Leonard and 33adianslcy for material which posseses two 
different wave velocities but his inputs were such that 
there was no discontinuiby along the slovrer wave fronts, 
fillers |4l applied the method of characteristics to solve 
a semi -infinite cylindrical (Timoshenko type) shell subo ac- 
ted to dynamic loading at one end. His input was also such that 
there was no discontinuity along the slower wavefront. 
Subsequently Chou and Koenig \ 3 \ apolied the raethod of 
characteristics for several spherical and cylindrical dila-ta- 
tion wave problems. Recently Chou and Mortimer l6| prooosed an 


* Numbers in parentheses represent reference in Appendix III 



unified approach- for one dimensional wave propagation 
problem wliere propagation of discontinuity along waves may 
exist. He proposed a method of solution for a system of n 
second order partial differential eq_uations with generalised 
displacement as dependent variables. 

1.3 OBJSOl AID SCOPE OF 'THE STUDY 

It IS felt that a thorough knowledge or the transient 
stresses developed in a dam-like structure subjected to 
ground motion will be useful for bhe advancement of theore- 
tical kno¥/ledge about the behaviour of such structure 
during an earthq.uake. For this study information from the 
records of a recent earthc[uake (11th December, 1967) at 
Koyna, situated at the south-west end of the Deccan trap, 
have been used. The physical details of the dam at Koyna are 
used in the numerical example of this study, .Analysis for 
'che longitudinal component of horizontal strong motion is 
Ignored since the dam is strong enough in longitudinal 
direction. To date, very few records of actual strong motion 
are available in this belt, and there is therefore little 
statistical basis for the prediction of further eq_rthquakes. 
Hence, the result of the analysis predicts the transienb 
stresses developed in Koyna-concret e-dan during earthquake 
on 11th December, 1967. 



Tlie problem is formulated, in Giiapter II, liie 
governing equations considering the hydrodynamic efCect are 
developed and characteristic equations as well as compati- 
bility along the characteristics are established. 

Recurrence relations for numerical step by step 
integration are given in Chapter III, Stability and conver- 
gence of numerical integration scheme are discussed. 

Thirdly in Chapter IV, summary of the investigation is 
outlined and conclusions are drawn, A comparitive study of 
the effect of change of cross section of a cantilever beam 
for a particular input ground motion has also been conducted 
in this section. It is felt that analysis developed herein 
will help visualising the transient stresses developed in 
structure subjected to ground motion. The analysis presented 
herein will be useful for the analysis of tall stack-like 
structure under dynamic loading. 

Detsmls of the matrices, determino.nts and vectors 
introduced in Chapter II and Chapter III are given in 
Appendix I. 

Finally in Appendix II, propagation of discontinui- 
ties in the domain of solution are discussed. The equation 
of oump across wavefronts are developed. Details of the 
problem solved by Leonard and Budianshy jlj are also given 


here 



CHAPTER I' 


POmiTriA-TION OF PROBLEM 

2.1 BISIC liQiJATlOiTSjIPITTAI APB BOUliDARY OOTOPTIOHS 

Concret'3 gravity dans c^re generally cast in. mono- 
liths. Tlie goint; between two subsequent monoliths are pa,cked 
with waterproof maberaals to prevent waterleahage. For 
reduction of ilermal sbresses in the structure it is custo- 
mary to avoid 1 .y interconnec bion of two adjacent monoliths 
(Figure l)t After bhe eai’chquahe on 11th December, 1967 
independenb de Coma' icn of aO^ecent monoliths of Koyna- 
concrete dara I'^’d been noticed clearly. Hence, bhe dam can 
be assumed for analysis as a cantilever beam with distribu- 
ted mass, according to the cross section of dam, along the 
height of the dam for a particular monolith. Since monolith 
has considerable widtb in transverse direction compared to 
the height of dam, shear deformation cannot be neglected, 
■Station of the cross s'^ction of dam will be small. However, 
for accuracy of anolys- '! correction for robation has been 
considered. 

Genera] assumotjcn for the material of the dam and 
impounding nut or arc the fcllowungs: 


(i) The mab'^'^ial of Lie dam (concrete) is linearly elastic 
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with VISCOUS daiup ing proportional to a constant/^ 
which varies, for concrete, from 0.1 to 0.15 times 
the fundamental critical damping 0. Fundamental 
critical damping corresponds to the limiting case 
of preiodic oscillation involving deformation of all 
the masses toward the same side and cixi he written as 

c c 

Fundamental critical damping C = 2 A 

( 2 . 1 ) 

where, = density of concrete 

0 

A = area of concrebe at a particulem cross 
section of monolith, 

C- ' ^ = fundamenbal freq_uency of the system. 

(ii) Effect of crack in the concrete is neglected, 

(ill) The impounding water is incompressible with irrotatio 
nal flow of small amplitude. The internal viscosity 
and some ether causes of energy dissipiation from the 
impounding water are disregarded, a.s well as the 
effecc of surface wave, 

2.1.1 Eg uat ion pf __Fib rat lo n o_f ^^B eam 

Assuiiing relative deflection w from bhe support 
and rotation Tjr of the cross section to be small bhe dynamic 
equilibrium ol Timoshenko beam (Figure 2), sub^jected to 



s 

ground acceleration Z^^(?igure 3), including the eilect 
01 the inertia and the shear deformation may he 

given as the followings: 


Id + t:i xjj = 0 


( 2 . 2 ) 


V - kA°G(w^ - t) = 0 (2.3) 

jL 

"‘x P°'^ tt - <■ = 0 (2.4) 

h - + 2tt) - « n = ° (2.5) 

v/here 

G ^ /U 0 

J = rotary inertia about y-axis. 

Other symbols are used in staxidard form. 


In the set of Eiuatims (2.2 to 2.5) subscripts 
refer to the p.^rtial differentiation with respect to 
variables mentioned there. 


For a cantilever beam subjected to su(3den ground 
acceleration the initial condition can be written as: 


w(x 50 ) - w^(x,0) - ■[|■(x,0) 


0 ( 2 . 6 ) 



and bne boundary conditions are 

M(L,t) = V(L,t) = w(0,t) 

= l|r(0,t) = 0 (2.7) 

Systeui of Equations (2,2 to 2.5) represents 
dynaroic equilibrium of a mjnelith of dam subjected to an 
eaarthquake where reservoir basin has no inipcunding v/ater 
v/ith initial conditions (Equation 2.6) ano boundary condi- 
tions (Equation 2.7), 


2,1,2 G-ov_erni n<_ E qu at ion o f Monolith of fail Vibration 

of the R serv o ir Bcsjui 

Hydrodynamic pres._ure will develop on the up- 
s bream face of a dam during 3.a earthquake in opposite 
phase with ground acceleration. Hence, equivalent inertia 
force of bhu mass of moving v/ater vvith the dam is to be 
added with bhe inertia force 'jf dai itself. The width of 
water b viiig with the dam during an earthquake in 
accordance with ganger's j7! theory is given by 


b = 




2 


( 1 - --w) + ( 1 

H 




( 2 . 8 ) 


where 

X - height of particular cru ss-section from bne 
base of the dam 

H = height of impounding water at upstream of 


the dam 



0^ = pressure coelfrcienb depending on the angle 

of inclination of the upstream face v\rith 
vertical. 

Hexice, mass of the vvater along per foot vyidbh of che dam 
IS given by, 

mass/it = A = : b (2,9) 

since unib width i/i the longitudinal direction is considered 
in this analysis. In Ec[uation (2.9) 


= density o f water 
= area of water 

Neglecting the hydrodynamic effect due to deforma- 
tion of dai, the governing differential equation of dyna.mic 
equilibrium of the monolith of a dam may be given oy system 
of Equations (2,2) to (2.4) and 



^tt “ ^^t " ^ ^ ^tt "" ° (2.10) 

where 

(-’A = 

The system of Equa.tions (2.2 to 2,4) ana Eauation 
(2.10) are the governing equations with initial conditions 
(Eauation 2.6) and boundary conditions (E ‘nation 2,7) for 



tiie vibration of a monolitii of dam subjected bo an earthquake. 


2.2 DERIVATION OF GOJIPATIBILITY CONDITION 

The input acceleration for response la bhe strucbure 
due to an earthquake is a seismograph reading and it is 
difficult to obtain any closed form representation of the 
saj'ie. Hence numerical solution of the problem will be sought 
in this sbudy. To facilitate numerical inbegration in conti- 
nuous formulation, it is advantageous to cha-nge the system 
of Equations (2.2) to (2.4) and(2.10) along with initial 
conditions (2,S) aad boundary conditions (2.7) in a physical 
plane (x,t-plane) to the characteristic form. 

Differentiation of equation (2.2) and (2,3) with 
resp ct to b yields 


+ El ijr^ = 0 (2.11) 

= 0 ( 2 . 12 ) 

EquatiOiis(2. 11 ), (2.12), (2.4) and (2,10) form a. 
system cC four simultaneous partial differential equations 
which in the non-dimensional form are given as: 


M-^ + I a- = 

V'c - 
Yi ■" ^ kj 


0 


3 / — 


(2.13) 


0 


(2.14) 




liie nondimensional pai'ameter a = 1,0 when there is 
no water in reservoir basin or when hydro dynojuic eilect of 
water is neglected and a = 1,0+ f A / p A for vibration 
of monolith o [ dam with a reservoir basin of water. 


The initial conditions given by fignat ion (2,6) change to 



v/(x,0) = v{^,0) = if(x,0) 

= a(-:E,C) = 0 (2.17) 

and the bjandar^;^ conditions (Euaabvjn 2.7) in nondimen- 
sional fjrni are 

M(i. r ) = v(i, r ) - -(0, r ) 

= t1i(0. D = 0 (2.18) 

A characteristic diiectioi of the system of Actua- 
tions (2,13 to 2.16) allov/ possible discontinuity oj the 
first derivative of dependent variables (a, v, iu, V) along 
its normal direction. Given in_tial data on such a direc- 
tion the normal derivative of the variables are not 
uniquely determined by the system of Equations (2,13 to 
2.16). Hov/ever, for legions in n, V --plane where the 
dependent variables are continrous the following relation 
must be satisfied: 


dM = 

M- 

dx 

-h 

M ^ 

d r 

(2.19) 

da = 

% 

d- 

X 

+ 

ap- 

dr 

(2,20) 

dv = 

% 

dx 

+ 

v^ 

d r 

(2.21) 

dV = 


dx 



dr 

(2.22) 



Equations (2.13 to 2.16 and 2,19 to 2.22) forai a 
s^stera of ei'di'c equations which, may be used to solve for 
eight first derivatives of dependent variables if the dis- 
tribution of dependent variables are knov'm along an initial 
curve. The system of Equations (2,13 to 2,16 and 2.19 to 2,22) 
can be given in matrix form : 

M [z; = [ff] (2.23) 

vdiere 

LJ -*■ represents a matrix 
[ I -♦ represents a vector 

The matrix X, vectors and Z are given in Appendix I, 

Sol''^iiig for M-^ from Equation (2,23) by well known 
Cramer’s rule yields 

= -4— (2.24) 

C |X| 

^/here 

1 I represents a determinant, 

jyj IS formed by replacing., one column of 
according to Cramer’s Pnle. 

Since first derivative of dependent variables may 
be discontinuous, the derivative is undetermined, A 

necessary and sufficient condition for the indeteruiinacy of 



IS [xj “ 0 and jx] = 0. 'I'iie vanisiiin'; of 

|xl leads aftea? applying Laplaoe ©acpansion teoimique nf 

determinants bo: 


l(dr )^ - J(dx)^ 





1 0 A A 


(2,25) 

'liierefore, eibJaer 


= t (2.26) 

dx ' 


dt Wo 

= + (,A/r= )V2 = + Op (2,27) 

dx ■■ ■ 


Hence the original s'/stem of Equations (2,2 to 
2.4 and 2.10) are totally hyperbolic since the roots of 
Equation (4.25) are real and distinct. The roots + 

(i = 1,2) are called physical characteristic direction and 
it IS custofliery to call inverse of it as the wave velocity. 

The Ccnisliing of |y| yields, the ch.- racteristic 
equation v/here d7( /dx is replaced ty + 0 | 


7dt ~ J da + dM = 0 


(2.28) 


Ihe equation (2.28) will be called coaipatibilicy 
relation along + and-G^ respectively, dt can be very 



or a- v/ill 


easily shown solving for eitiier M~ or a-^ 
result in sets of equation whicli ai e idencical to Liqua- 
tion (2.26) a,nd Equation (2.28). 


Similarly solvin,-? for 7-?- or V- or or v~ 

will lead CO the identic=.l equation-compabibility relation 
along + Gg. liter replacing df /die Toy + Og the compati- 
bility relation along + are given as : 

- a + + |3v! 1^"+ '} dv — ^ = 0 (2.29) 

~ ^ I 


where 

i 

> 



A 




V 



Solution of the dependent variables may be oi^banaed 
in the region, v/here the dependent variables are continuous, 
by solving the system of four Equations (2.28 to 2,29) 
simultaneously along the characteristic directions (Lfua.- 
tion 2.26 o.nd 2,27) with proper initial conditions (Equation 
2.17) and boundary conditions (Equation 2,18), 



CHiPT^R III 


FJj.II RIG aL SOLiEIOF 


3.1 DERIIAIIOTJ OR EaCJ {.RBi'CE RSLA'L'IOI 

Solu’cion of tiie syoteiu of equabions (2.26 to 2.29) 
cax'i be oboai.aed niinierically bf seep by sbep inbegration 
vfaicli amouiibs to maa'ching out of solution, from 'TT = 0 

by means of Hartree [sj technique of numerical intS'^ra-tion. 
In perfonning numerical integratiorij tiie fixiicc length in 
the X direction is divided by strips at an interval h 
(J ~ 1, J, J + 1 etc. in Figure 4) v^hich inturn form a 
mesli network with the constant time line (R ~ 1, K, K + 1 
etc. in Figure 4). Ihe dependenb variables (a, v, V, M), 
in this process, are considered as known functions of x 
at a time , either as given initial conditions or as 
the output of previous stage computations, 

let the dependenb variables bo knoi/m on a line 
Tr" = and is required to be found on a line X.~V^ +21 cT 
on which a typical pivotal point is J, K (Figure 4). 'There 
c.ce four laailies of characteristics (+*^^ and +O 2 ) passing 
ti-irough e3,ch interior pivotal point (grid). These families 
of characteristics incersect the initial line K-1 at four 
points (I-l, I~2, 1+1, and 1+2 in Figure 4 ) ~ which hence- 
lorbh will oe called subgrids. Since coordinates of pivotal 







point JsK IS laiown, using liquation (2.26) and ±jquation 
(2.27) spocs coordin-.te x of subgrids may be obtained, 

Eo'ijGion (2.26) in finite diCference forn. reduces 

to : 

- &^7o^ (3.1) 

Since is the slope of the characteristic line 
and for a rectangular cross section (l = J) it has a cons- 
tant slope. 

Similarly, Eq_uation (2.27) yields 



^1-1 

= X. - 4t 

J 

/02 

(3.3) 

a.nd 

^1+1 

= 

/Cg 

(3.4-) 


In equations (5.1 to 3.4) subscripts refer to the 
subgrids OT grids. 

Since the dependent variables (a, v, ll, Y) axe known 
functions at tliree ad]acent mom grid points (J, J-1 and J+1 
111 Eigure 4) on the initial line 7-1, knowing the coordinates 
of m subgrid a quadratic interpolation formula as given 
below' may be used to find each of the variables. 



_ _ X - _ 


(- _ j )2 


2li 


2” 


(3.5) 


<±LGrc a , ^ . refers bo the functional value of a 

3+ I , K- I 

at frid point J+1 and on the line K--1 and x refers to the 

J 

s.jace coordin?i.te of j. Ihe nomencla-ture introduced will he 
used hereafber in this study. 


Fox.'iulfie simil.j.r to Equation (3.5) ca.n he ’.uabten 
easily for the remaining dependent variables. 


Eoi a p. rticul r pivotal point J, K the couipat ihility 
relations ( uq.uations2. 28 and 2.29) are written in finite 
difference for i in terms ob the knov/a values of dependent 
variables at suhgrids. 


Ecpiation (2.28) re -laced hf finite difference gives 


along +0^ : 


7 ^^3,k V2,k-1^ 


- ■2- Ji-2^®'3,k - ®-i-2,k-1^ 


- V2,k-l) 


(3.6) 


and along 


0 



^ V2,lc-1 


)dr - ^ " ®'i+2,k-1 


+ " V2,lc-1 


) = 0 


vjiicre 


'^1+2 ' ^2 ^ 


and ^x-2 '^D ^ ^x-2 


C5.7) 


Siiiilp'nly E mat ion (2.29) reduc-s 


CO 


along +C2- 


_^)dC + ( A + t:^) T 


D,lt 


^ n-i.k 

( a' - 

= _ ^(a + oc^_^)d'^ (3.8) 


and along - Cgi 
^ /r 


P d 




^dr 


V, 


+ ( > - ^^)^i+1,k-1 ■" " ''1+1, k -1 


) 


= F(a + a^_^^)dt' 


(3.9) 


Ax-1 = 


where 



and 


'^ 1+1 ^ -^3 Vl 




(fk.1 + f^)/4 


The refers to groand oc cole rat ion at time t!" = K 

ixid a refers ro value of the sam’i ec grid or subgrid men- 
D 

tioned in subscript. 


In v.;rii;ing Equations ( 3. 6 to 3.9) it is assumed that 
Zl’Zr IS suiricientlv Siiail so th-'b linear relation holds 
between iho pi'/otal pome and initial line ^ nd hence y avera- 
ging of depun j-oiit variables may be done. 

The system oe Equations (3.6 to 3.9) may be written 
in matriEE uot-^cion as bho following: 


\;heLU 


[b] (u) = [c] [pj + j_Qi; (4.10) 


|b] a 4 id jo] = 

= 




CoePficient matrix 

vector of the dependent varic-bles 

at pivolal point J,K. 

vector representing dependent veri- 

ables at subgrids 

vi-ctor of input ground acceleration. 


Equation (3.10) may be reduced to give a recurrence 


r..! non for step by step integration for pivotal mesh point 



cis encouxito ued, whicli i& ^iven by 


{u; = [sr'' Icl [pj + [b]''' {q} (3.11) 

where 

[bJ””' = inverse of LbJ . 

Befalls of all vectors and matrices used in liqua- 
tion (3.11) aj e given in A.ppendix I. 

For bliG ve.tical boundary plane (x = 0 and x = 1) 
suitable recurrence relation may b- obtained using proper 
boundary conititions and following the same general princi- 
ple OL inteiior mesh point. At the left boundary where 
X = 0 (Figure 5) two families of characteristics (- 0.j and 
- Co) intersect the bnundarr pivotal point where the 

boundary conditions from Equation (2.18) ai^e 

1(0, r ) = v(o, r ) = 0 

I'ho space coordina.te subgrids ( 1+1 and 1+2) may 
be obtainc-d from Emation (3.2) and Equation (3.4). From 
l-cnown values of dependent variables at initial line, the 
functional value of variables at any subgrid m may be 
obtained by using an interpolation formula as the following 



7 






- v2 

(x , . - X ) 

+ .ai , (V 


2ii 


3+2, k~1 “ ^^^3+1, k~1 




(2.12) 


Formulae sx'nilar bo Enuation (5.12) may be used for 
teraaining dependent variables. Introduction of knov/n boundary 
coxiditions in finite difference form of compatibility relation 
along - 0^ and - C 2 yields: 


along - 0^ 




Y 


1+2, k-1 


__^ ) d^ + i J 


7 '^ 1+2 ^- 1 + 2 , k-1 


^1 “ \+2,k-1^ ° 


(3.13) 


and along - G 2 ‘* 


^1+1, k-1 + ( ? - ^■^r^)'^i+1,k-1 




(3.H) 


A recurrence relatiOxi for left boundary may be obtained 
solving Enuation (3.13) and Aquation (3.14) simultaneously 
and given by. 



I] = [d] ^ [H] [r| + [nr'' (s? (3.15) 


where 




= vector of dependent variables as speci- 
fied at pivotal point J,iv 


iR| = vector of dependent variables at subgrids 


= vecbor of the input ground acceleration 


fs] 

Lt and I R'j are coefficient matrices. 

Details of matrices aiicl voc cores used in Equation 
( 3 . 15 ) are given in Appendix I. 

Sijnil.<-rly for right boundary x = 1 (Figure 6), 

where boundary conditions from (2,18) are 


M( 1 , r ) = V( 1 , T ) = 0 

Ihe spcLce coordinate of subgrids are found from 
Eauation (3.1) and Equation ( 3. 3 ). Dependent variables at 

■4“ lo 

m subgrid ijioy be evaluated from the quadratic interpola- 
tion formula as given in the following ; 


X 


X 




" %k-1 ~ ^3-2, k-1^ 


D-1 


- ^2 

. ^^ 1-1 - 


la. 


+• a 


D-2,k-P 


(3.16) 



Pormulcie simil-ir to E5_uation (3.16) laay "be used 
lor other dcpendenb variables, 

A recurrence relation sinilar to Equation (3,15) 
rarr/- be obtained hy substituting proper boundary conditions 
in Equations (3,6)and Equation (3,8) and is ^,iven b/: 


whore 


and 






K (5.17) 


y~j and |n] 


UJ 


Djk 



coefficient raat rices 

dependent variable vector at pivotal 
point J,j<; 

vector of dependent variables at 
sub grids 

vector of ground acceleration. 


Details of matrices .^nd vector used in Equation 
(3.17) are '..ivoen in Appendix I. 


Tie solution oi the set of Equations (2,26 to 2,29) 
n 0.7 be obtained using the recurrence relations (Equations 
3,10 , 0.15 and 3.17) end proper initial conditions 

( Iquation 2,17). 



3-2 STABILITY AID COI/GflfiriGE 


To illustrate the general uniqueness theorem |9 1 of 
linear hyperbolic system of equations, let be a non- 
characteris 'cic line (Figure 7) of continuous slope along 
vhiich the dependent variables are prescribed as continuous 
functions. A unique solution to the system of equations 
(2,2 to 2,4 axia 2,10) may be obtained in bhe triangular 
region/ p cut out by the outer characteristics (+ 0^) at 
an interval in T from and having an intercept on 

it. The interval on the line is bhe domain of 

dependence of P on Hence, it is implied that the ratio 

o-r step size ( / Zi 5) should comply with the outer 

characteristics. 

The stability of finite difference scheme with rectan- 
giiloa: grids \va,s studied by Courant, Preidrichs and Dewy |lOj 
for a wave equation (a™ = u^^ ) and the condition of sta- 
bility and convergence towards the analytical solution as 
A X and approach zero was found to be 

^ 1 (3,18) 

A X 

In the case of wave equation, characteristics are 
tvvo lines (d71/dx = domain of depen- 

dence of P IS /^p (Figure 8). So, iT a step size is chosen 





sacii that AT/ Ax > 1, the pivotal point V then lies 

beyond the doiiiain of dependence ^ reached by" bhe continu- 
ous solution. Hence, such a, solution will hardly converge to 
xhe analytical one. However, accuracy of the solution is 
adentical with bhe analytical one v^hen AT/ Ax = 1 
since in tnis case pivotad. point coincides with P, Accuracy 
of solution vfill decrease with decreasing ratio of step 
size owing to f'Ocrease of slope of the lines 'VP^ and YS 2 
which play the role of finite difference cha-racteristics . 
However, a physical interpretation can also be given very 
e isily. If for instance, the step size is taicen less Than 
bhe inverse oC bhe faster wave velocity, it only amounts to 
artificial f£i,ster propagation of th? wave. Hence, a step 
size is to be chosen as to comply with faster wave velocity 
for convergence oC solution to The analytical one, 

Ihc proof of convergence of bhe Hartre ^ scheme wa.s 
also given b" Sauer [nj xn 1952. 

5.5 HTJiAKiGiA IHTlCTR/fflOH 

The height of the da.a u divided by main grids at 
equcH inTcrval. ( dx = h). The step size in /T direction is 
always chosen a.ccording to JSnaation (2.26), During integra- 
tion, it IS assumed that bhe accelerogram reading is linear 
between Oxie peak to subsequent peak. The time inteiwal 
between subsequent peaks is divided in n equal parts of 



step size h. a,xiJ. remaining portion At; loriis another 
step, He/icej n+1 stej:is are required dor integration of 
one peak io other, The procedure for step hy step inte- 
gration 'Vhen dt = dx = h can be summarised as follows: 

St ep ( 1 ) ; 

The space coordinate of subgrids are located using 
Ec[uations ( 5. 1 to 3.4). 

St ep ( 2 ) : 

Since = 1 for the problem, the dependent vari- 
ables at interior subgrids (Figure 4) ^^1 and 1-1 are 
interpolated from the known variables at main grids on 
initial line using suitable interpolation formula similar: 
to Eouation (3.5) or Equation (3.12) or Equation (3.15). 

Step (3): 

De|iendent variables at main grids are then calcula- 
ted using suitable recurrence relations from Equation (3.10, 
5.15 and 3.17). 

For the last step AC( step (4) and step (2) 
remain identical except step (2) whece variables at all sub- 
grids have to be interpolated. 

By continuou^y repeating the basic computational 
steps the solution may be marched out from ^ = 0. However, 
for starting integration at ^ = 0 initial conditions 
(Equation 2.17) have to be used. 



OHAPTJiR I? 


DISCUS Slop JLTD COBCLUSIOF 

4,1 DISGU33IOU OIT iUJ145RICAL A-PAIYSIS 

Stabilitiy and conTergence of tiie numerical sciieme 
proposed in bins SGudy has been discussed in Chapter III, 
I he \ral 1 d 1 t 7 ol the seheme is shown by coinparing the 
resalbs obtained for a iroble'U with a Icnown closed— form 
sola bion. The conrergence of the numerical scheme is also 
seen by obtaining solution of this problem for various 
grid sizes. 

The numerical scheme has been applied bo obbaiii 
the vari tioi of base shear force and base noment with 
bime (Figure 9) of a cantilever Tiaaoshe )l £0 Beam vnbh 
molerial property such that tho flexural and shear waves 
ha.ve ssuTt® velocity (H 0 ) ax'id subjected, to um lines/r 
velocity ( V ) at the base. The numerical resulbs are in 
good agreemenb with that obtained by Leonard -ind Sudian— 
shy hi by Laplace transform technique, llcror for bhe 
tv/enty-one grid sizes is not more than one unit. This 
pcoblem has been given in detail in Appendix II. 


Subsepuejibly, a tapered cantilever oea-i (0 - 1^) 
hs been sub 3 Gcted to a step input fer diffirenb durations 



of 1:11116 of tiie ground acceleration, Hie corresponding varia- 
tion of base shear force and bace moment v^iitii time are 
evaluabed vo tli the help o£ the numerical technique developed 
in Chapter III (Figure "lO and 11 ), Peaks of bhe response 
?/ith time sho'/ diff-rences for various grid sizes for bhe 
input of short duration tuc the general trend of the curves 
are quite simiLar, On the o bher hand, when the dura,tion of 
ground accelerobion is large^ a little deviation in bhe 
shear-time ojj-'rc is obnerv 3 d at network of twenty-one 
grids from th'*' shear-birae curves oC finer networks (Jl and 
81 grids), flic moment -time curves are same for all practical 
purposes (diff:rence in fourth significanb figure). 

In aay numerical a^nalysis the accuracy of the solution 
is D_imited bp truncation and round off errors. Furtliormore, 
a meaningful solution vail be obtained if and only if the 
solution converges as :hc grid size is reduced. In bhe case 
of small duration of ei.c taaon the disparity in pea„k vatLues 
for various grid~sizos may be due bo the fact that the 
response crea/ted on the structure by bhe ground acceleration 
IS apprecicbLy small, S^nce the step by step integration 
uses the previous stag“c values of the variablos to predi.ct 
the same at an inmerval of time, there is possibility of 
error accumulation v/lr'u bhe previous stage results are not 
quite accurate, I/[oreo\cr, fur coarser grid aizes bhe 



possitHity oi incurring large error cannot be ignored aiid 
bins ej^laiiis the flattened response, Hov/eiyer, when bhe 
duration of exitation is large the solutions show reasonable 
convergence for diffirent grid sizes (Figure 11). Sifnilarity 
of results obtained by Leonard and Budiansky jl | and that 
obtained by nuxierical approach establishes perhaps conclu- 
sively the stability and convergence criteria of the 
solution as well as accuracy of the sarne, 

4-2 DlSObSSXOii and GOlTGLUSIOh 

Firstly^ change in the response of cantilever beams, 
with varying cross-section, subjected to ground acceleration 
has been studied. The material of the beam is assumed to be 
linearly elascic with no damping. The reduction of base shear 
force and base moment occur as the angle of taper (0) 
increases (Fiyu^^e 12). This follows directly from the 
elementary theory because reduction of the mass of the 
beam must give less moment and less shear at bhe base. 

Secondly, effect of damping has been considered for 
3 . camtilevsr beam of linearly elastic material. The variation 
oj. base shear force and base moment with time has been 
plotted li'i Figure 13 and 14. The curves show that the damp- 
ing ratio has a predominant effect on the response of 
structure us the time interval increases. The wave length 
of bhe response also increases with increase of the deuaping 



rabio. ' IS in good agreement with, the theory " 

A reduction of two tnird la the peak -values 
of the bending moment and shear force at base are obtained 
betweer the cases when the daii^ing ratio is zero and when 
lb IS 0,1, In this analysis the absolute damping is assumed 
proportional bo fundamental model freq_uency due to inherent 
lifliitation of knowledge on the absolute damping of the 
material. This may lead to a prediction of higher responses 
from actual in the transient state of solution. However 
any firm conclusion on the prediction of actual response 
c uonot be given unless the absolute daraping of material is 
known with reosonable accuracy. 

The method is finally employed to predict the 
transient stresses developed in the highest mcnolith of 
the Koyna-Concrete-Dam during the earthquake of 11th 
December, 1967. The material of the dan has been assumed 
linearly elastic with two diffirent viscous danpiii-; rstios 
( = 0.0035 and = 0.1). ihe absoluce damping of 

the material lias been evaluated using the fundamental 
frequency as given in an unpublished report with the help 
of numerical integration approach. The variation of bending 
stresses and shear force with time are plotted in Figure 15 
to 38 at tliiee different cross-section of the dam (A, B and 
0 as shown in Figure 1), The following two cases are analysed 
in this study; 



(a) The response oJ the dai having no T/ater in the 
reservoir basin (Figure 15-26). 

(b) The respoxise of the dara with water in the reservoir 
basin (Figure 27-58). 

In the former case (a) ihe inaxiniuiii transient response 
reduces about Cifty per cent as tne damping ratio increases 
from = 0,0035 to - 0.1. The maximuji transient 
response increases al'iiost by sixty per cent in the latter 
case (b) from that of the former (a). In both the cases 
rapid changes in the values are indicated from the results, 

'Then reservoir basin has no water and damping ratio 
OJ the material of the dam is 0.0035 the bending stresses 
crceed the allowable tensile limit at all cross-sections of 
the dam. These results have only theoretical importance of 
estimatiiig the error thit might occur in the response if 
the material ol the dam is assuxiied bo be linearly elastic 
without VISCOUS dajrpingc Hovvever, when the damping ratio is 
0.1 the predicted maximum bending stresses exceeds the 
allov/able liu-'ib only at section *G' which coifor’ji to the 
cracking o bh. sectioa at that level. 

The maximum bending stresses are more than the 
allowable stresses in all the three cross-section o the dam 
for both values of damping ratio considering full wacer in 



tiae reservoir bctsin., Th. stress predicted are qui'ce high 
owing to biie li.^iitaticn of Zanger's \l\ theory. The concept 
of addition oC virtual nass is liable to predict maximum 
stresses on che higher side than that of actual because 
the mass of water moving with the dam will not remain 
constant foe the v^-hole interval of time. Moreover no energy 
dissip ration ii the water is considered in this stud3r. 

In a strong earthquake the bendin stresses may be 
su CCicientll y large to erceed. allowa.ble stresses for a brief 
period oi cimc. This does not imply failure oC the structure 
since the sc cesses will be restrict d as soon as the ground 
motion and/or hydrodynamic force reverse or decrease. 

However, it ^.s felt thac further study should be 
directed tov^/ards the generalisation of ma,terial properties 
for the da i. Since, concrete is neither elastic nor elnsto- 
plastic and the stresses abeam a ver/ high value only for 
a short duration of time, investigation should be carried 
out Gonsiderni- the time dependent material properties and 
nonreversib] e factors such as redistribution of stresses 
due to cracking, energy dissipation due to daraping etc. 

Aga,in due to enormous bulk of a dam, the accurate predic- 
tion of matecial behaviour v;ithin it is quite impossible. 
Hence, p co ^sbilist ic approach may be very suitable for cliis 
type of study. 



Sxnce tvro dimensional behaviour Ol sub^ecbed 
bo irouiid accelerabion r^nd ixiteraction beb?/een vibrating dam 


luid the impounding water as a composite s'/sbem, 
thoroughly e nlored, investi.'sations may also be 
towards these lijoes. 


has nob been 
direc bed 













FIG. i2 EFFECT OF CHANGE OF CROSS SECTION ON 
VARIATION OF SHEAR AND MOMENT WITH TIME 
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iiPefiuDix - I 


AI,1 D-JiXiiLi OF MaTRICXS, DFI0R4I , j'li' j ^ . j \7 lJC1’ORS HjJD 
Il'T Ai'.ul"SIS 

The de bails o-' che biatricesj det ninaats a.id 
YOCbors used lu Chapter II and Ghapcer III a:3 giveri iii 
the follOYinas: 
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dt 3x 0 0 0 
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dr dx 0 
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J t dx 
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1 -> 0 
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(A1,3) 
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appendix - II 


BO'DID^^RI C0NDIII0F3 INYOIiVUiG B'£ SB IiIpTJI 

a2.1 PROPaMIOF of DI soon in TIT I 

liic: solution of the sustain of E Tua.t ions (2.2- to 
2.-4 and 2,10) vvith bounda-ry conditions Involying step 
change in dependent variables (a, v, ?, S) presents 
complications of the type discu-sssd in GJnapter II, For 
oth. -r boundiiry conditions imposed eJ-ong each of the two 
lateral bound,, ries of t,he corr-3po.nding serni-infinite 
strip in the .half plaAe >/ 0, t.he solution to t.b.e above 
on a-n x~inte:cval #1 finite length does not picesent si-iilar 
problems. The normal derivatives of dsp-anclent variables 
across a physical characteristic may be discontinuous^ 
but t.hese di-scont inuities do not affect the application of 
Aqiiations (2,26 .and 2.21). However discontinuities in 
dependent variables are not governed bv Equations (2,28 
and 2.29). This necessitates the deter aina-t ion, firstly 
of t.he lino in the half space bounded by finite interval 
ox X, along w.hich the discontinuities occur, and seconbAy 
the magnitude of such discontinuities. 

That the discontinuities of dependent vocciofoles 
occur along characteristics may be proved by lassaming; .a 



non-chaj?act eristic line / J across which cilscontiaait 7 
e;?:ists i.e., b-jtwean and (figure It is 

furcther as3'a:.a .d that in such a rbgioa, bounded hf 
a.ncl aJ..i. functions are continuous, Integratio of 

d uation (2.29) witn. upper sign along +O 2 from point 1 
to 3 gives 

ry 3 

f PdZ + A (v, - vS - / — - = 0 (A2.1) 

1 ^ ^ 1 I 

where 

* ^ _ 11 ^ 

Q “ - a +- -T \iv 

Assuming Q is bounded, as 7^2 tends to p^ i.e. d'Z. -O, 
dquoubion (A2.1) reduces to: 


= 0 (A2.2) 

A 

Where S stanuo for the jump of .a pajrbicular vaxiable 
specified. 


Similarly, integration along -Gg -s -» 0 gives 


/\6v + 


M = 0 
I 


(il2.3) 


Gomtaination of Equations (A2.2 and a 2.5) yields: 


6v = 


6? 


0 


(A2.4) 



ilong +G^ as 6 


0, it ca,n 


I'ollowing siaiilar procedure 
be sho'TO. b.ha/b 


6M =* 6a = 0 (A2«5) 

The iD ;uations (a2.4 ana A2.b) show that the ^innp 
across a noncharacteristic line cannot exist, 

A2.2 GOhDITIOr OP JUMPS 

The condition of jump across a chsiract eristic e.g 
may be obtain .^d by writing E -ucicion (2.28) with lovver 
sign for points 1 and 2 (Figure 12-2), The following 
Eciuation showing condition of jump along +C^ is obtained 
point 1 tends to point 2 

Along +0 ^ : 



Siiailarly conditions of jujup along the other ciiaracter- 
istics may be given as: 

ilLong -C^ : 

6M = - J (A2.7) 


Along +02’* 



6Y 


(A2.8) 




JO-ong -02^ 
6? 


( A^"- « 9 ) 


k 6V 

The Jii -■nations (a2,6 to A2.9) reveal that firstly, 
discontinait 7 of M is associated 'Hha’': that of a aaid dirs™ 
continuity of Y is associ a+.ed with t. hat of v, and secoiicilys 
discontinuities of M and V are unco'u'pled. 

A2.3 MA&IIfUlh 0? JUMP 

'T,h.e vniignitude of gum'p across tC,| ma3r be obtahiierl 
by writing .Scraation (2.28) with u-upsr sign for two lines 
and /~2 (with one b-low and one a,bove +G,j) onid. 
subtracting one from other, l-hus, as /^.| tends to 
gives 


- J d(6a) - 0^ d(6M) = 0 

Elimination of 6a by E-uation (a2.6) yields 

j- d(C. 6 M/J) + d(6M) = 0 

Oi 1 

which after differentiation reduces to 

_ d( 6M) 

6M 


d(j/C^) 

~/G^ 


1 

7 


(A2, 10) 



Solution of ji]niation (A2.10) may be given by 




i[Z 

a/Oi^j 


(A2.11) 


Wliere (J/g ^)2 indicates value of the same at point 2 
wiiich. is the upper limit across +C^ and the corresponding 
lovver limit is 1 across +C^. 


In a similar way the magnitade of jump across 
from point 1 to 2 may be obtained* ‘I'his is given by 


V 


2 


1/2 


V- 



(A2, 12) 


It com be very easily verified tho.t the jump 
across -G^ will be identical with bhe Equation (A2,11) 
a,nd jump across --C 2 will be identical as the Equation 
(A2.12). 


A2.4 UflEORa BEid! fllH = tig 

leoimrd o;nd Budiansky jlj solved a pro bleu of 
uniform boom subjected to a jump input v(0,^ ) = 1,0; 
t >/ 0, It v/ces .assumed that the rant e rial of beam i.. 


linearly elastic and possesses a property such that both 



waves have the same velocities of propagation (kG- = Ji) 


The ch;araot eristics epuation from -triU'.tions (2,26 and 
2,27) are as the following: 

“ = +1 = (h2.13) 

dx ~ ~ 


and the condition of jump along the char-acteristtcs from 
Equations (a 2.6 to 42.9) a,r8 : 


Along +C^ : 

5M « 6 a (.42.14) 

67 = ~ > 6v (42.15) 

Along -C^ : 

6® = - 6a (42.16) 


67 = j^br (42,17) 

The .magnitude of j-ump across +G^ reduces to 
Along +C.J : 

1 =^ (6S) - 1/2 (67) = 0 (.12.18) 

(6V) + /■/2(6a)= 0 (42.19) 



Solution of the system of .iouations (A2.18 ond Aii.19) are 

6a = A* sin ( 7’ •“ c' ) (A2,i0) 

67 =» A* QOS Z - (,12.21) 

Where A and G are constants of integration which are 
to be determined from the boundary conditions as prescribed 
in the problem,. I'or the particular problem, for which the 
boundary conditions are ^(O/C ) 1 and ) ~ 0 for 

T yo, the constant ax integration hav^ the following 
voAues 

A* = - 'Tf 

and G =0 

Hence, Bqu-otions (A2.20 and ii2.2l) reduce to 

sin > S/2) (A2.22) 

and 6? =* - / cos ("'^7' x/2) (A2,23) 

Thorefare jump of the dependant v ria-bles odong +0,^ is 
given by 



6a(x, x) 


= .rr 


y sin x/2) 


(12,24) 



6V(x,x) 


> Cos ( '' Z' x/2) 


(A2,25) 


The boundary conditions prescribed ab the other 
end are T ) " V(1, '^ ) =0 for V 0. Gornbination 

of _-icpiatioas (a 2.14 and a 2»24) yield 

61(1,1) = 6S(1,1) = -r>“BinrF^/2 

6 ?( 1 , 1 ) = cos f “^"“/2 

Hence, to satisfy the boundary conditions a. ^nriip 
hf varia.bles v/ill •ccur along -G^ such that 


6a(x, 2-x) = 

- 6M(x, 2-x) = 

- sinifpT x/2) 

(A2.26) 

6?(x, £-x) = 

7- 5v (x,2-x) = 

7.cos(V^ s/2) 


(A2,27) 


In a 

siiail-.j,r 'w 

ay it can be 

shown that 



5 7(x, 

2+x) = 

^C93(V7' 5 

;/2) 

(A2. 

,28) 

6a(x, 

2+x) = 

sin('f7 

x/2) 

( A2 . 

, 29 ) 

6V(x, 

/i--x) = 

- 7 . cos(Y^ 

5/2) 

(A2. 

, 30 ) 

6a(x, 

4--X) = 

V 7 sin(Y> 

x/2) 

(A2, 

>31) 



Beyond = 4, it can be very e.isily verified tlia-t 
each intar-7^,‘.l of =4 will have simalar jumps identical, 
to dpuat ions ( .k 2 ,24 to A2 , 3 1 ) . 

A2,5 NlMeRIOvLL ItOHNIQUfi 

ITuinerical techni-iue developed in Chapter III nia-y 
also be used for the types of problem involving ste:) chanpe 
of dependent v riables. Ho-.-ever, afditioneJ- co-nplications 
A^ill arise duo to tracing of 7ia,ve fronts (Giioa'^acteriscics) 
and addition o.f proper jumps along the vvave fronts. 
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